# TODO: Add comment
# 
# Author: guochun
###############################################################################
library(spatstat)

sigma2=alpha=numeric()
for(i in 1:10){
	X <- rLGCP("matern", function(x, y){ 1 - 0.1 * x +0.1*y},
			c(mean=0.0, variance=2, nugget=0.0, scale=3, nu=2),
			win = owin(c(0, 10), c(0, 10)))
	re=kppm(X, ~x+y, "LGCP", covmodel=list(model="matern", nu=2))
	sigma2[i]=re$mcfit$par[1]
	alpha[i]=re$mcfit$par[2]
}
